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ABSTRACT 

We study how uncertainties in the rate coefficients of chemical reactions in 
the RATE 06 database affect abundances and column densities of key molecules 
in protoplanetary disks. We randomly varied the gas-phase reaction rates within 
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their uncertainty limits and calculated the time-dependent abundances and col- 
umn densities using a gas-grain chemical model and a flaring steady-state disk 
model. We find that key species can be separated into two distinct groups ac- 
cording to the sensitivity of their column densities to the rate uncertainties. The 
first group includes CO, C+, H+ H2O, NH3, N2H+, and HCNH+. For these 
species, the column densities are not very sensitive to the rate uncertainties but 
the abundances in specific regions are. The second group includes CS, CO2, 
HCO+, H2CO, C2H, CN, HCN, HNC and other, more complex species, for which 
high abundances and abundance uncertainties co-exist in the same disk region, 
leading to larger scatters in the column densities. However, even for complex 
and heavy molecules, the dispersion in their column densities is not more than 
a factor of ~ 4. We perform a sensitivity analysis of the computed abundances 
to rate uncertainties and identify those reactions with the most problematic rate 
coefficients. We conclude that the rate coefficients of about a hundred of chemi- 
cal reactions need to be determined more accurately in order to greatly improve 
the reliability of modern astrochemical models. This improvement should be an 
ultimate goal for future laboratory studies and theoretical investigations. 

Subject headings: accretion disks - astrochemistry — methods: statistical — 
molecular processes — stars: planetary systems: protoplanetary disks 



Introduction 



More than 140 organic and inorganic molecular species consisting of up to 13 atoms have 
been identified in s pace so fafl including such con aplex species as di methyl ether (CH3OCH3; 
Snyder et al.lll974l ) and acetamide (CH3CONH2; iHoUis et al.ll2006l ). The rich variety of the 
observed molecules implies that many more relevant yet undiscovered species must be in- 
volved in the processes of their formation and destruction. Still, only a tiny fraction of these 
species have been firmly detected in several protoplanetary disks, including C O (and its iso- 
topes), CN, HCN, HNC, H2CO, C^H, CS, HCO+, H ^ 3qq+^ DCO+, and NoH+ JPutrev et al. 



19971 : iKastner et al.l 119971 : Ivan Zadelhoff et al.ll200ll : lOi et al.ll2003l : iThi et al.ll2004h . Multi- 
molecule, multi-transition observations of emission lines with (sub-)millimeter interferom- 
eters and single-dish antennas provi de a wealth of info r mation about planet-forming disks 



Petv et al. 


2006; 


Dutrev et al. 


2007) 



^http : //astrochemistry . net| 
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An ultimate goal of disk studies is the reconstruction of the evolutionary history and spa- 



(Willacv & Laneerl 


2000; 


Markwick et a 


2002 


van Zadelhoff et al. 


I2OO3 


Semenov et al. 


2005; 


Aikawa & Nomura 


2006; 


Tscharnuter & Gaill 


2007) 



20031 ; llkner et all 12004 



20071 ). The analysis and 



modeling of observational data acquired with limited spatial and spectral resolution is not 
complete without taking all possible uncertainties into account such as instrumental errors, 
uncertainties in the distance to source, orientation, etc. Whereas these uncertainties can in 
general be reduced either by using better observational data provided by current interfer- 
ometers (PdBI, SMA) equipped with a new generation of receivers or in the future with the 
advent of more powerful instruments like ALMA, EVLA, and Herschel, there is an intrinsic 
source of ambiguity in all the chemical models - our limited knowledge of the reaction rate 
coefficients. 

Modern astrochemical da t abases include up to about 4500 gas-phase reactions and y 
450 species jMillar et al.lll997l : lie Teuff et aPboOol : ISmith et al.lbooi IWoodall et aDl2007h . 
but many of these reactions have poorly estimated rate coefficients with uncertainties of 
about a factor of 2 and larger. For example, radiative association reactions lead to the 
formation of new complex species out of smaller ones through the photon relaxation of an 
excited coUisional molecular complex. Under low-density space conditions where three-body 
processes are unlikely, these reactions may have rate coefficients as high as ~ 10~^ cm^ s~^ 
and as low as ~ 10~^^ cm'^ s~^, depending on the density of vibrational states of the complex. 
The radiative association rates are usually difficult to measure in the laborato r y as well as t o 
predict theoretically, especially for bigger species (lBateslll95ll ; IWilliamslll972l ; lHerbstlll980l ). 



In contrast, dissociative recombination of molecular ions almost always proceeds very 
rapidly, in particular at low temperature, with rate coefficients that can be accurately ob- 
tained (IFlorescu-Mitchell fc Mitchellll2006l ). However, the products and their relative branch- 
ing fractions branching channels of many dissociative recombination reactions are not known 



precisely, particular ly at low temperatures (see though ISemaniak et al.ll200ll;llshii et al. 



OiekuU et al.l 



Geppert et al. 



20061). which may spoil theoretical predictions (see, e.g. 



Millar et al. 



2005af ) 



2006'; 
19881 : 



Recently, IVasyunin et al.l (l2004j . Paper I hereafter) investigated the influence of un- 
certainties in rate coefficients on molecular abundances in dense and diffuse clouds, using 
the gas-phase UMIST 95 database. They found that the abundance uncertainties of simple 
species are limited to within about an order of magnitude and increase substantially with 
the number of atoms in the molecule, though the uncertainties still do not exceed the obser- 
vational errors for simple molecules. They proposed a sensitivity analysis t o identify those 
reactions that introduce the largest errors in the computed concentrations. IWakelam et al. 
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( 120051 . hereafter Paper II) studied uncertainties in the gas-phase chemistry of hot cores and 
demonstrated that at late times, > 10^ yr, modeled abundances of impo rtant molecules 



can suf fer from large uncertainties due to poorly known rate coefficients. IWakelam et al. 



(l2006al ) focused on dense cloud chemistry and took uncertainties in the physical parameters 
into account, as well as comparing the osu.2003 and UMIST99 chemistries. These authors 
came to conclusions that grain-surface reactions should be considered in order to achieve 
better consistency with the observational data, even though there is a large degree of am- 
biguity in simp l e cloud models based on static physical and chemical structures. Finally, 
Wakelam et al.l (l2006bl ) expanded this study by using their uncertainty analysis to show 
that there is a bistability in the abundances of many species that are hyper- sensitive to the 
ado pted value of the ratio of the cosmic ray ionization rate of helium to that of hydrogen (se e 



also Pineau des Forets et a. 



1992; IShalabiea fc Greenberd Il995l : iBoger fc Sternberg! 120061 ). 



Furthermore, llzzard et al.l (l2007l ) recently studied the effect of proton-capture reaction rate 



uncertainties in the NeNa and MgAl chains on the abundances of the Ne, Na, Mg and Al 
isotopes produced in intermediate-mass AGB stars. 

Despite previous work in rate uncertainties, to the best of our knowledge, there have 
been no attempts in the literature to study how these reaction uncertainties, which are 
unlikely to be eliminated in the near future, affect the results of astrochemical modeling 
in the wide range of physical co nditions typical of pro toplanetary di sks. In this paper, we 



extend the previous analyses by IVasyunin et al.l (120041 . Paper I) and IWakelam et al.l ( 12005 



ipi 





Paper II) to the conditions of a low-mass protoplanetary disk surrounding a young T Tauri 
star. The aims of our study are several. First, we analyze how large the abundance and 
column density scatter for key species are due to uncertainties in rate coefficients. Second, we 
investigate how these uncertainties vary with disk location. Third, we isolate the reactions 
whose rate uncertainties contribute most to the abundance scatter and which are therefore 
worth studying in detail. Last, we predict how the overall consistency of theoretical models 
will be improved after these reaction rates are better constrained. 

The organization of our paper is the following. The adopted disk physical structure 
and time- dependent gas-grain chemical model are presented in Sect. 12.11 The Monte Carlo 
method used to introduce the uncertainties in the RATE 06 reaction rates is described in 
Sect. 12.21 In Section [3] we report computed distributions of the mean molecular abundances 
and column densities and their errors in the whole disk, and analyze the influence of the rate 
coefficient uncertainties on these quantities. A correlation method that allows determination 
of the role of various reactions in the abundance uncertainties as a function of time and disk 
location is outlined, and results for several key species are presented in Sect. 14.11 Particular 
attention is paid to the identification of the most uncertain reactions for the chemical evolu- 
tion of key species in the entire disk. In Sect. 14.21 we predict to what extent the abundance 
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scatter due to the rate uncertainties will be lowered when accurate rate values of several 
dozens most problematic reactions are available. The remainder of Section H] is concerned 
with problems in different classes of reactions. Final conclusions are drawn in Sect O 



2. Disk model and uncertainty approach 
2.1. Physical structure and chemical model of the disk 



In our simulations, we adopted the 1+lD steady-state irradiated disk model with vertical 
temperature gradient that repres ents the low-ma.ss Cla ss II protoplanetary disk surrounding 
the young T Tauri star DM Tau (jP'Alessio et al.lll999l ). The disk has a radius of 800 AU, an 
accretion rate M = 10~^ M^r^ y^".^; ^ vi scosity parameter a = 0.01, and a mass M ^ 0.07 Mq 
(jPutrey et al.l 119971 : iPietu et al.l 120071 ). The thermal and density structure of the disk is 
shown in Fi g . [H A s a disk age we used a value of ~ 5 Myr, which has been derived by 



Simon et al.l (120001 ) based on evolutionary track modeling of the central star. In chemical 
simulations an outer disk region beyond the distance of 50 AU from the central star is 
considered. This is the only routinely accessible disk region with existing (sub-)millimeter 
interferometers. 

We assumed that the disk is illuminated by UV radiation from the central star with 
an intensity G = 410 Go at r = 100 AU (G (r) oc r~°'^) and by interstellar UV radiation 



with i ntensity On in plan e-parallel geometry (jPraind 1 19781 : Ivan Dishoeckl Il988l : iBergin et al. 



20031 : iDutrey et al.l 120071 ). The dust grains are assumed to b e uniform 0.1 /xm sph erical 



particles made of amorphous silicates with olivine stoichiometry (jSemenov et al.ll2003l ). with 
a dust-to-gas mass ratio of 1%. The self-and mutual-shielding of CO and H2 against UV 
photodiss o ciatio n is computed using the pre-calculated factors from Tables 10 and 11 in 
Lee et all Jl996h . 



Three other high-energy sources in the model that drive chemistry either through dis- 
sociation, ionization, or desorption are cosmic rays, decay of short-lived radionuclides, and 
stellar X-ray s. The X-ray ionization rate in a given disk region is computed according to 
the results of iGlassgold et al. Jl997al jbl) with parameters for their high-metal depletion case. 
In this model, the thermal ~ 5 keV X-ray photons are generated at ~ 0.1 AU a bove the 



central star, with a total X-ray luminosity of ~ 10^^ erg cm ^s ^ (IGlassgold et al.ll2005l ). 



The cosmic-ray ionization rate is assumed to be 1.3x10 s ^. Th e ionization rate cause d 
by radionuclide decay (primarily ^^Al and ^"Fe) is 6.5 x 10^^^ s~^ (IFinocchi fc Gaill 119971 ). 



We adding the X-ray and radionuclide ionization rates to those involving cosmic ray particles 
(CRP) and CRP-induced UV photons. 



- 6 - 



A modified static gas-grain chemical model that inclu des gas-grain interacti ons without 
surface reactions and turbulent mixing has been utilized ( Semenov et al.l 2005). We allow 



only s urface formation of molecular hydrogen using the approach of iHoUenbach &: McKee 
( 119791 ) ■ This modification significantly decreases the numerical demands of our chemical 
code and makes it computationally tractable while allows us to include all necessary physical 
processes. The gas-phase reaction rates and their uncertainties are taken from the recent 
RATE 06 database, i n which the effects of dipole-enhanced ion-neutral rates are included 
( IWoodall et al.l 120071). To cal c ulate photoreaction rates through the disk, we adopt pre- 
computed fits of Ivan Dishoeckl (119881 ) ins tead of integrating the wa velength-dependent cross 
sections over the local UV spectrum (see Ivan Zadelhoff et al.ll2003l ). 



Gas-grain interactions include the accretion of neutral molecules onto dust surfaces with 
a sticking efficiency of 100%, dissociative recombination of ions on charged dust grains, 
and grain re-charging, as well as UV-, CRP-induced and thermal desorption of surface 
species. Desorption en ergies E^es are mostly taken from the recent osu.2007[^ database 
(iGarrod fc Herbstll2006l ) or roughly estimated by analogy for about 20 molecules. 



Overall, our network consists of positively and negatively as well as neutral dust grains, 
420 gas-phase and 157 surface species made of 13 elements, and 5773 reactions. Among 
these 5773 reactions are 4517 gas-phase reactions, 2 charge exchange reactions for dust 
grains, 940 dissociative recombination reactions of ions on charged and neutral grains, and 
314 accretion/desorption processes. 



As initial abundances, we have adopted the so-called "low metal" set of iLee et al.l (119981 ) 
as listed in Table [1] and assumed that all hydrogen is mostly locked in its molecular form. 
In these initial abundances, the standard solar elemental composition is depleted in heavy 
elements by, e.g., 200 times for S up to more than 10^ for Fe, CI, P, and F. The mostly 
atomic initial abundances are chosen instead of those from a molecular cloud model because 
this choice allows the abundance uncertainties to accumulate in a long sequence of chemical 
reactions with imprecise rates, s tarting from basic pr ocesses that reach steady-state at early 
times during the disk evolution. IWillacy et al.l (119981 ) have shown that most of the molecular 
abundances are hardly affected by the choice of input abundances. Note also that the 
resulting abundance un certainties can be sensitive to the initial elemental concentrations 
(IWakelam et al.ll2006bl ). However, we do not allow the elemental concentrations to vary in 
order to keep the computations in reasonable limits. 



^https : //w ww . physics . ohio- state . edu/~eric/reseaxch . html 
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2.2. Method to model rate uncertainties 



Each reaction rate coefficient in the RATE 06 database is given in the standard Arrhe- 
nius form and thus relies on 3 parameters: a (the absolute value at room temperature), 
/9 (the index for the power law dependence of the rate on temperature), and 7 (the acti- 
vation energy barrier in K). A specific expression for the reaction rate coefficient depends 
on the type of chemical process: bimolecular reaction, direct cosmic-ray ionization reaction, 
photodissociation, etc. Note that in all cases the rate coefficient scales linearly with the 
parameter a. 

In addition to these three parameters, the RATE 06 reaction rates are characterized by 
an accuracy estimate (A, B, C, and D), where the uncertainties are smaller than 25%, 50%, 
within a factor of 2, and within an order of magnitude. Almost all reactions with measured 
rates belong to the ffist group (~ 1 400 reactions), while those reactions with rates that were 
"guessed" or derived by analogy are assigned to the third group (~ 2 800 reactions), which 
contains most of the ion-neutral reaction rates. The second group includes about 300 mostly 
neutral-neutral and photodi ssociation reactions, whereas the fourth group consists of only 4 
photodissociation reactions (IWoodall et al.l 120071 ). 



In accord with the latest osu.2007 network, we adopted a rate uncertainty of one order 
of magnitude for radiative association reactions. Such a large uncertaint y factor is ju stified 



Herbstlll985h . The 



by the difficulty in calculating or measuring these rate coefficients (e.g. 
uncertainty factors vary widely for the photoreaction rates since the corresponding frequency- 
dependent cross-sections have limited accuracy and have been obtained only for a fraction 



of molecules in the UMIST database (Ivan DishoecklllQSSi : Ivan Dishoeck et al.ll2006l ). 



Our method to model abundance uncertainties is based on the computation of a large 
set of chemical models, using identical physical and initial conditions and the same chemical 
network, but with randomly varied rate coefficients within their uncertainty limits. We 
utilized the same method to introduce uncertainties in the rate coefficient values as described 



m 



Dobriievic fc ParisotI (ll998f ). lDobriievic et all (l2003f ). and lWakelam et all (l2005f ). The rate 



coefficient for each gas-phase reaction i is randomly chosen over a log-normal distribution 
with median and dispersion Fj. Consequently, the rate coefficient of the i-th reaction 
spans an interval between ttj/Fj and aj x Fj with a probability of 68% (lo"). We generated 
a sequence of these log-normally distributed rates for all gas-phase reactions in our chemical 
network using the following expression: 



where a- is the /-th realization of the z-th reaction rate, is the standard RATE 06 rate 
value for the i-th reaction, Fi is the dispersion, or the uncertainty factor, for this rate, and 
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e is randomly distributed from -1 to 1 with an uniform distribution law. 

Using this approach, we allowed the rates of only 4517 gas-phase reactions to vary in our 
network. The rates of all other processes (gas-grain interactions, dissociative recombination 
of molecular ions on grains, and recharging of dust grains) were kept constant. Not did we 
account for possible variations among the parameters of the disk physical model. Such an 
idealization allowed us to focus solely on an investigation of how uncertainties in gas-phase 
reactions affect the results of the disk chemical modeling. 

Using 8 000 realizations of the RATE 06 network, we simulated 5 Myr of chemical evo- 
lution in the outer, r > 50 AU, disk of DM Tau. The disk grid consists of 5 radial (50, 100, 
200, 380, and 760 AU) and 10 equidistant vertical points (with step sizes of 3.2, 7.6, 18, 
41.7, and 96 AU for the considered radii). Such a huge number of utilized chemical models 
with randomly varied rates leads to low values of statistical noise in the distributions of the 
abundance scatters at any particular disk location, assuring the correctness of the sensitivity 
analysis. Moreover, a large number of varied reaction rates allows an analysis of all reactions 
in the che mical network, in cont rast to models with limited chemistry included in previous 



works (e.g. JWakelam et al.ll2005l ). Still, even with the modest 5 x 10 spatial resolution of the 
adopted disk model and without surface reactions, the overall computational time needed to 
calculate the chemistry with these 8 000 networks was about 3 days on a 4-CPU PC machine 
(3.0 GHz Xeon, 8 Gb RAM). 



3. Results 

3.1. Abundance distribution profiles 

Before we perform a detailed analysis of the influence of the rate uncertainties on the 
modeled molecular abundances, we investigate the abundance distribution profiles. Because 
we utilized a log-normal distribution for the rate coefficients, one might expect that for each 
molecule the abundance distribution should show a normal (Gaussian) profile in logarithmic 
scale. If this hypothesis is correct, one can use the dispersion - a measure of the scatter 
of randomly varied values that have a Gaussian distribution - as a convenient quantity 
characterizing the abundance and column density uncertainties. 

As shown in Paper II, this assumption may not be fullfiled for certain species and 
certain time steps. For example, when molecular abundances of a particular species show 
a steep decline or increase with time then at that moment the corresponding histogram of 
the abundance distribution can have several peaks and thus be far from a Gaussian shape 
because of sampling effects. Bimodal distributions can also be obtained if the system is 
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highly sensitive to a small variation of specific parameter. IWakelam et al.l (l2006bl ) showed 
that dense cloud chemistry hypersensitive, and even bistable, to the ratio between He and 
H2 cosmic-ray ionization rates in particular conditions. In this case our uncertainty method 
would tend to overestimate the dispersion of the modeled abundances. The disk chemistry 
is very rich in a sense that it proceeds vastly differently in various disk regions and at 
different times, so it is natural to expect that the abundance distributions of some species 
can sometime deviate from a normal distribution. 

We carefully studied this problem and found that such a situation happens rarely at 
late times, > 1 Myr, for which our analysis will be performed. In Fig. [2] we show several 
representative histograms of the abundance distributions at 5 Myr, using HCO"'" and CO as 
an example. For an inner disk region, at r ^ 97.5 AU, 2; ^ 30 AU, a strong deviation from 
the Gaussian shape is clearly visible, with 2 distinct peaks and a gap in between (Fig. [21 top 
panels). Note that this is one of the few exceptional cases that is discussed in detail below. 
Everywhere else through the disk the HCO"'" and CO abundances have nearly Gaussian-like 
profiles (Fig. [2l bottom panels). This condition holds for other key species as well, implying 
that the averaged abundances in the disk should not differ from the values computed with the 
standard RATE 06 network. Moreover, in Paper I we found that the Gaussian shape of the 
abundance distribution is also preserved for the uniform distribution of the rate uncertainties 
in linear scale. 

In the rest of this paper, we will refer to one standard deviation (Icr) in the mod- 
eled abundance distribution as the "abundance uncertainty". In our notation "abundance 
scatter" represents twice the abundance uncertainty in logarithmic scale (2 log cr). It is the 
abundance uncertainties that have to be taken into account when interpreting the observa- 
tional data with a chemical model. 



3.2. Distribution of the mean abundances in the disk 



Although the evolution of all species was calculated in 8 000 realizations of the chemical 
model, in what follows we focus on several key molecules that are used as tracers of disk 
physical structure and chemical composition. Our hst of important species includes molecules 
that are widely used for kinematic studies and determination of the density and temperature 
(CO, CS, H2CO, NH3, and HCN) as well as for probing such major parameters as ionization 
degree (HCO"*"), ra diation fields (CN, HCN), and the yet-to - be detected C"*" as a dominant ion 



in th e disk surface (jPutrey et al.lll997l : iKastner et al.lll997l : lAikawa et al.ll2003l : iDutrey et al. 
2007h . 
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We use time- dependent abundances to analyze the chemical evolution of these species 
in detail. The analysis is based on the calculation of the relative importance of each reaction 
in the chem ical network for th e abundance change of a molecule under investigation at a 
given time (IWiebe et al.l l2003l ) . This allows us to weight all reactions in the model and 
isolate a few most important pathways for chemical evolution of one or several species under 
specific physical conditions. In the case of protoplanetary disks, such an analysis has to be 
performed for several tens of representative cells. 

In Fig. [3] we show the mean abundances of these molecules in the disk at 5 Myr, which 
are obtained by averaging the results of the chemical network. The vertical extension of the 
disk is scaled by one hydrostatic scale height, if (r) = \/2C,JVL, wh.eYe Cg is the midplane 
sound speed and VL is the Keplerian frequency (jPartois et al.ll2003l ). The unity pressure scale 
height is about 25 and 320 AU at distances of 100 and 800 AU, respectively. 

The "sandwich" -like chemical structure of the disk is clearly seen in Fig. [31 In the 
"dark", cold, and dense midplane at < 0.5 — 1 pressure scale height, the molecules freeze 
out onto dust grains within 100 — 100 000 yr and low steady-state fractional abundances of 
~ 10^^^ — lO^^^are reached. The abundances of the gas-phase species in the midplane at late 
times are sustained by the CRP-induced desorption of the mantle materials. The midplane 
abundance of CO ~ 10~^, is higher because this molecule is very abundant (~ 10~^) in the 
gas phase ini tially, and since CO i ce is rather volatile because of its low CO desorption barrier 
of ~ 930 K (IBisschop et al.ll2006l ). Overall, the gas-phase abundances in the midplane tend 
to slightly increase with radius as the disk surface density and thus the density in midplane 
decrease and gas-grain interactions become less intense, which slows down the freeze-out of 
gas-phase molecules onto dust grain surfaces. 

The surface layer of the disk {Z/H{r) > 2) is deficient in molecules due to strong irra- 
diation by high-energy stellar and interstellar photons. Even CO, which is self-shielded and 
mutually-shielded by H2 against UV photodissociation cannot survive. Among the consid- 
ered species, only C"*" reaches its maximal concentration and becomes the most abundant ion 
in the surface layer, with an abundance of ~ 10~^. The thickness of this highly ionized layer 
(~ 1 pressure scale height down from the disk surface) stays nearly constant with radius for 
all key species but highly reactive radicals CS and CN, whose abundances are already high 
(~ 10~^) at scale pressure heights as high as Z/H{r) ^ 2 a.t all radii (Fig. [3]). 

In contrast, polyatomic species tend to concentrate in the warm, chemically rich inter- 
mediate layer (< 1 — 2 pressure scale height, depending on the radius). The model-averaged 
abundances for the considered species, ~ 10^^° — 10""^, are close to t he values computed with 



the standard RATE 06 network and those f ound in previous work (lAikawa &: HerbstI Il999 



Willacv fc Langei]l2000l : lAikawa et al.ll2002l : IWillacv et al.ll2006l : iDutrev et al.ll2007f ). 
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Now let us consider the chemical processes that form and destroy the considered species 
in the intermediate layer. In contrast to dense cloud models, for many species in this disk 
zone steady-state is reached by ~ 10^ — 10^ yr. This is due to the higher densities encountered 
in the warm layer (> 10^ — 10^ cm~^) than in the cloud cores (~ 10^ — 10^ cm~^) and a non- 
negligible flux of the UV radiation. Imp inging interstellar UV photons allow many surfac e 



species to desorb back into the gas phase (IHasegawa fc Herbstlll993l : iGarrod fc Herbstll2006l ). 



thus reducing their surface abundances and overall importance of gas-grain interactions for 
disk chemical evolution at late times, t > 10 000 yr. 

The evolution of CO and HCO"^ in the intermediate layer is governed by a small set of 
reactions at all radii. The CO molecules are formed via reactions of atomic oxygen with CH, 
CH2, and C2 at t < 1 000 yr, and partly through the dissociative recombination of HCO"'" 
afterwards (at this stage steady-state has been reached and there is a loop between CO 
and its protonated form). The light hydrocarbons that are precursors of CO are themselves 
rapidly created within a few years via the radiative association of C with H and H2, and 
CH with H2. The CO is mainly removed from the gas phase by accretion onto the dust 
grain surfaces in outer, cold regions and by reactive collisions with helium ions in the inner. 
X-ray ionized part of the intermediate layer. In the upper part of this layer, CO is also 
photodissociated. The ion HCO"*" initially forms via reactions of O with CH and ionized 
light hydrocarbons, CH^ (n=l-3), and the ion- neutral reaction between CO and CH^ (also 
producing methane). At a later time of > 1000 yr, HCO"^ becomes the most abundant 
ion in the intermediate layer and the dominant molecular ion in the disk, with a fractional 
abundance of ~ 10^^. At this steady-state stage, HCO"^ is mainly produced via reactive 
collisions between CO and Hj]" and destroyed by dissociative recombination with electrons 
and negatively charged grains. 

The chemically related H2CO molecule is produced through reaction of CH3 with oxy- 
gen atoms. At a later time of ~ 10^ yr H2CO is involved in a simple formation-destruction 
cycle. It starts with protonation of formaldehyde by HCO^, HaO"*", and H^, which is fol- 
lowed by dissociative recombination into H2CO (33%), CO (33%), or HCO (33%). Apart 
from photodissociation and freeze out, another direct destruction channel for H2CO in the 
upper, more ionized part of the intermediate layer is the ion-neutral reaction with C^, which 
produces either CO and CH^, HCO+ and CH, or H2CO"'" and C. The chemical evolution of 
formaldehyde typically reaches steady-state at about 10^ yr with a gas-phase abundance of 
< 3 X 10^^°. Note that in the same disk chemical model as discussed above but with a set 



of surface reactions taken from lGarrod fc HerbstI (120061 ) the overall (gas-phase and surface) 



abundance of formaldehyde is decreased by a factor of several as it is converted to methanol 
by hydrogenation on grain surfaces. 
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The photostable radical CN is abundant (10~^) at an upper, less molecularly rich disk 
layer {Z/H{r) ~ 1.8, Fig. [3]). The initial formation channel for the cyanogen radical is 
the radiative association between C and N, followed by the neutral-neutral reaction of N 
and CH, and dissociative recombination of HCNH"*" at later times [t > 10^ yr). Other, less 
important formation pathways for CN at this stage are neutral-neutral reactions between C 
and NO producing CN and O as well as between C and OCN leading to CN and CO. The 
destruction of CN mostly proceeds through photodissociation and neutral-neutral reaction 
with N forming highly photostable nitrogen molecules and C. 

In turn, HCNH+ is produced by ion-neutral reactions between H2 and HCN"*" and be- 
tween N and C^H^ {n = 3, 4) at early times {t < 10^ yr). Later the steady-state recycling of 
HCNH"*" is reached, which involves dissociative recombination of HCNH"*" into either HCN 
(33%), HNC (33%), or CN (33%), followed by re-production of HCNH+ via protonation of 
HCN and HNC by either H^, HsO"*", or HCO+ as well as via the ion- neutral reaction of CN 
with H3". 

The chemistry of HCN is tightly related to the evolution of CN and HCNH+. HCN is 
mainly produced by neutral-neutral reactions involving nitrogen and methylene or methyl 
at early evolutionary stages, t < 10^ yr. Later the evolution of HCN proceeds as a part of 
the HCNH"*" formation-destruction loop. The major destruction routes for gas-phase HCN 
include freeze-out onto dust grains and the ion- neutral reaction with C"*", which forms CNC"*" 
and H. Another important destruction channel for HCN in the upper disk layers at ~ 1.5 — 2 
scale heights is photodissociation. Steady-state for HCN is typically reached within > 10"^ yr, 
with an abundance of about a few xlO"^''. 

The layer where abundant ammonia exists is located somewhat deeper toward the disk 
midplane compared with HCN and es pecially CN because NH3 can be photodissociated by 



UV photons shortward of ~ 1950 A( van Dishoeckl Il988l ). The main formation route for 



ammonia is dissociative recombination of NH^, which itself is formed through a sequence of 
hydrogen insertion reactions starting with the production of the NH"*" ions from N+ reacting 
with H2. Upon formation, NH^ repeatedly reacts with molecular hydrogen and gains an 
addition hydrogen atom until NHj is formed. The destruction of ammonia proceeds mainly 
via ion-neutral reactions with C"*" and S"*" in the upper intermediate layer as well as Hjj", 
HaO"*", and HCO+ in a deeper region, which leads again to NH4 and either H2, H2O, or CO. 
For all radii, the abundance of ammonia reaches a steady-state value of < 3 x 10~^° — 10~^ 
at very late times of > 10^ — 10^ yr. 

The chemical evolution of CS is governed by neutral-neutral reactions between sulphur 
atoms and either CH, CH2, or C2 initially, and later, at t > 10^ yr, by an extended formation- 
destruction cycle of HjiCS"*" (n=l-3) and HC2S+. At this evolutionary stage these complex 
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ions are produced by ion-neutral reactions of protonated water, Hg , and HCO+ with HnCS 
(n=0-2). The initial routes to form the HCS+ ion are the reactions between CH^ and S or 
CS^ and H2, whereas HC2S"'' forms predominantly via reaction of CH3 and S+. The H3CS+ 
ion is mostly produced in reactive collisions of methane with sulphur ions. Upon formation, 
HCS"*" disso datively recombines into CS (~ 10%) or CH (~ 90%), while HC2S~'' breaks up 
equally into C2S or CS. The H2CS+ and H3CS+ ions recombine with equal probabilities into 
CS and either HCS or H2CS. Finally, the major destruction channels for CS are photodisso- 
ciation and the ion-neutral reaction with C"^. The steady-state abundance for CS, ~ 10^^, 
is usually attained at about 10^ — 10^ yr in our model. 



3.3. Distribution of the abundance and column density uncertainties 

For each key molecule, we derived the distribution of abundance uncertainties (disper- 
sion) in the disk. The results are shown in Fig. |H These uncertainties are quite different 
in different disk regions. Note that the overall abundance uncertainties are not much larger 
for the more complex (heavy) species among considered, in comparison to those with fewer 
atoms (e.g., compare C"*", CO, HCO"'', and H2CO). However, if we would focus on even more 
complex species (e.g., consisting of > 7 atoms), the abundance uncertainties wil l in ge neral 



increase with the numb er of atoms in the molecule, as shown by lVasyunin et al.l (120041 ) and 



Wakelam et all (!2006al ). 



Among the key reactions, the rates of photodissociation and photoionization processes 
can be uncertain by factors of 2 and larger. M oreover, their absolute valu es depend on the 



details of the UV radiative transfer modeling (Ivan Zadelhoff et al.l 120031 ). Thus, it is not 
surprising that in the disk atmosphere at Z/H{r) > 1.5 — 2 large abundance uncertainties 
with a factor of < 10 — 30 are reached for most of the species, in particular across the interface 
between the intermediate layer and the disk atmosphere. A notable exception is ionized 
carbon, because it contains essentially all of the element carbon in the disk atmosphere and 
thus maintains a high steady-state abundance that is not affected by chemistry. This effect 
leads to relatively low values of the C"*" abundance uncertainties, in particular a factor of < 3. 
The same is true for CO, which locks up all carbon in the disk midplane and intermediate 
layer at < 2 pressure scale height and here has its lowest abundance uncertainty, a factor of 
<2 (Fig.SD. 

There is a disk region adjacent to the midplane at r ~ 100 AU, Z/H{r) ~ 1, where 
the abundance uncertainties of the considered carbon-bearing species including CO show 
an increase of up to one order of magnitude (see HCO"*", Fig. H]). There, the chemistry of 
CO and hence all chemically related species (e.g., HCO+, H2CO, etc.) are prone to the 
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hypersensitivity at > 1 yr caused by the rate uncertainty of the X-ray ionization of hehum 



atoms, as discussed in IWakelam et aLl (j2006bl ). This region is moderately warm, with a 



temperature T ~ 25 — 40 K, and well shielded from the UV radiation from the central star 
and the interstellar UV radiation [Ay ^2 — 5 mag). The thermal bremsstrahlung X-ray 
photons (~ 5 keV) generated at ~ 0.1 AU above the star, however, are able to penetrate into 



this region (iGlassgold et al.lll997al ). with a total ionization rate of ~ 10~^^ s~^. These X-ray 
photons produce helium ions that rapidly react and destroy the CO molecules (forming O 
and C+). In turn, these C+ ions are used to reform CO but also to slowly produce carbon 
chain molecules (e.g., C2H2) and cyanopolyynes (e.g., HC3N), which are removed from the 
gas phase by freeze out onto the dust surfaces. These heavy surface molecules lock up most 
of the elemental carbon at t > 1 Myr such that the gas-phase abundance of CO is reduced 
to ~ 10~^ — 10^^. The overall efficiency of the carbon chain formation and the decline of the 
CO abundance at late times sensitively depend on the He ionization rate. Consequently, the 
CO histogram of the abundance distribution at 5 Myr has a double-peaked, non-Gaussian 
shape, as shown in Fig. [2l 

The abundances of nitrogen-bearing species such as CN, HCN, and NH3 in the region 
around the midplane at r > 200 AU (Fig. H]) possesses a high sensitivity to the rate of the 
CRP ionization of He. The histograms of CN, HCN, and NH3 are not fully Gaussian at 5 Myr 
in this region, and our method tends to overestimate the resulting abundance uncertainties. 
This region is so cold that at late times, t > 3 x 10^ yr, most molecules are removed from 
the gas phase by freeze-out. At this evolutionary stage, H;^ and later H+ become the most 
abundant charged species. The primal formation pathway for H+ is the slow ion-neutral 
reaction between He"*" and H2. The rise of the H"*" abundance at ~ 5 Myr(< 10^^) occurs 
in part because of the slowness of the radiative recombination of ionized hydrogen atoms 
in the molecular-deficient gas. The total electron concentration increases in tandem with 
H"*", and leads to more rapid dissociative recombination of HCNH"*" and NH4 (precursor ions 
for CN, HCN, and NH3). Consequently, the previously low gas-phase abundances of these 
nitrogen-bearing molecules increase at t ~ 5 Myr by factors of several. The steepness and 
evolutionary pattern for these abundance gradients depend sensitively on the abundance of 
helium ions. 

In the rest of the intermediate layer, the chemical evolution of the gas-phase species is 
initiated by radiative association reactions with high uncertainties but governed later by ion- 
neutral and neutral-neutral reactions with many accurately estimated rates. The situation 
in the disk midplane is similar, though at later times the chemical evolution of many gas- 
phase species is governed by the steady-state accretion-desorption life-cycle, which has no 
uncertainties in our model. This effect leads to lower abundance uncertainties of a factor of 
~ 3 than occur in the upper disk region (see Fig. H]). 
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It is interesting to note that for CO, C"*", and NH3 the uncertainty peak does not coincide 
with the maximal relative abundance in the disk. One might expect that the scatter in their 
column densities could be smaller than for other considered species (CS, CN, HCN, HCO^, 
and H2CO), for which uncertainty peaks and maximum abundances tend to coincide. 

In Fig. [5l the radial distributions of the column densities at 5 Myr are plotted for 
the considered species (dotted lines). In addition, in Table [2] we compile typical intrinsic 
uncertainties of the column densities of a larger set of key molecules in protoplanetary disks 
at 100 AU. Indeed, the overall uncertainty in the column densities of CO, C"*", and NH3 as 
well as , H2O, N2H"'", and HCNH"*" is less than a factor of about 2. Using this uncertainty 
value as a criterion, we assign these species to the so-called first uncertainty group. On the 
other hand, the column densities of CS, CN, HCN, HCO^, and H2CO as well as CO2, C2H, 
and HNC have larger error bars up to a factor of 4. These species belong to the second 
uncertainty group. It is likely that the species from the first sensitivity group are more 
reliable observational tracers than the species from the second one. 

Note that the column density uncertainties do not exceed a factor of about 4 even for 
formaldehyde and are comparable with observational uncertainties. However, for heavier 
species (e.g., HCgN), which are not included in our study, the uncertainties are likely larger, 
as shown in Paper I and Paper II. In our model, these uncertainties slowly increase with time 
and typically reach steady-state at > 10^ yr. The column density uncertainties can be used as 
error bars of the theoretically predicted quantities when comparing with observational data. 
Moreover, one can envisage a situation when, for an object with a well-studied structure, 
high-resolution observations followed by advanced modeling will allow putting such tight 
constraints on molecular abundances or column densities that it may help to determine the 
rate coefficients of some key reactions observationally. 



4. Discussion 

In the previous section, we demonstrated how and to what extent the uncertainties of the 
rate coefficients affect the computed abundances and column densities of several key species 
in the disk. In this section, we discuss the relative contributions of individual reactions to 
the uncertainties of these key molecules. 
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4.1. Sensitivity analysis method 

The chemical evolution of even a simple molecule usually involves a large number of 
reactions in the adopted chemical network. The non-linearity of the equations of chemical 
kinetics often makes it a challenge to find a direct link between the rate value of a certain 
reac tion and the m olecular abundances. Thus an efficient numerical method is required (see, 



e.g., lTuranyilll997l ). 



Using a sensitivity analysis based on a simple correlation method, we have shown in 
Paper I that under the physical conditions of dark and diffuse clouds the rate uncertainties 
of only a handful of chemical reactions affect the accuracy of the resulting abundances to the 
highest extent. For the disk chemistry studied in this paper a modified, two-stage sensitivity 
analysis was performed. 

First, for each selected molecule, we calculated time-dependent linear correlation coeffi- 
cients between the abundances and rate coefficients for all 8 000 realizations of the RATE 06 
network, and 19 logarithmically taken time steps over the 5 Myr of disk evolution at each 
disk cell. The coefficients Rij{t) are given by the expression 



R^Jir, z, t) = =, (2) 



Ez i^iir, z, t) - Xi{r, z, t)) J2i ("j - "i) 

where x\{r,z,t) is the molecular abundance of the i-th species computed with the l-th real- 
ization of the UMIST network, Xi{r,z,t) is the standard (mean) abundance of this species, 
aj is the j-th reaction rate coefficient for the l-th realization of the RATE 06 network, and 
is the standard rate coefficient of the j-th reaction. 

This is in contrast to Paper 1 where the analysis of statistically significant reactions was 
performed for only the final time step. Since at various evolutionary stages the chemistry of a 
molecule is typically governed by a restricted number of chemical processes, at any particular 
time step the abundance correlation will be high only for some of the key reactions, as shown 
in Fig. [5] for the case of ammonia. The obvious drawback of such a "single-time" approach 
is that it does not account for the cumulative character of the abundance uncertainties. 
Therefore, we integrated the absolute values of time-dependent correlation coefficients over 
the entire evolutionary time of 5 Myr. These quantities are called "cumulative correlation 
coefficients" . 

With this cumulative criterion, one can identify those reactions whose rate uncertainties 
are the most important for the abundance scatter of one molecule in one specific disk region. 
As a general criterion for the entire disk, we utilized the cumulative correlation summed up 
over all 50 disk grid cells, followed by summation over the 8 considered species (CO, C+, CS, 
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HCO"'", H2CO, CN, HCN, NH3). One should bear in mind that the results of our analysis 
rely on the inherent assumption that the absolute values of the rate coefficients and their 
uncertainties in RATE 06 are accurately determined. This assumption may not be fuUfiled 
for all important reactions in the RATE 06 database as we will show below. Such reactions 
will not be treated correctly by our sensitivity analysis. 



4.2. Identification of the most problematic reactions 

We calculated disk-averaged cumulative correlations for all 5773 reactions in our chem- 
ical network and for 8 key species: CO, C+, CS, HCO+, H2CO, CN, HCN, and NH3. The 
corresponding cumulative correlations, normalized in descending order, steadily decline with 
the number of reaction. The 56 most "problematic" reactions and their weights for these 8 
molecules are listed in Table [31 We term by "weight" the corresponding cumulative correla- 
tions summed over the all 50 disk cells and renormalized to the maximum value of 1148.1, 
which is achieved for the RA reaction between H2 and C. 

These 56 reactions correspond to 1% of the total number of reactions in our chemi- 
cal model. Moreover, such an apparently small number of reactions is chosen to stimulate 
experimental and theoretical studies of the rate coefficients of the reactions with high astro- 
chemical importance. Note also that some of the most "problematic" reactions in Table [3] 
have been carefully studied as their rates have small uncertainties of only factors of 1.25-1.5. 

In order to demonstrate that these 56 reactions are indeed important and their rate 
uncertainties strongly affect the abundances, we re-calculated the disk chemical evolution 
but allowed the rates of these 56 reactions to vary by only a factor of 1.25 (and 2 for radiative 
association). Such idealization is thought to mimic the situation when all of these 56 rates 
will be accurately measured in the laboratory or derived theoretically. 

The resulting uncertainties in the molecular column densities are shown in Fig. O (solid 
lines) and reported in Table [3] (last column). Compared with the initial column density 
uncertainties, the refined chemical model leads to smaller error bars for most of the considered 
molecules, in particular for CO, HCO^, H2CO, CN, and H2O. However, the improvement is 
not so obvious for nitrogen-bearing species (HCN, HNC, NH3, HCNH+) as well as CO2, CS, 
and C2H. Their chemical evolution is governed by a large set of key reactions that are only 
partly included in the list of the 56 most problematic reactions (Table 

The fact that the uncertainties in column densities decrease not only for the 8 key species. 



■^An extended version of this table is freely available from the authors upon request. 
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but also for some other molecules in Table [2] clearly indicates that some of the reactions from 
Table [3] are relevant for their chemical evolution. For these are the ionization of H2 and 
reactions of H3" with C, O, and CO. As we have shown in Sect. 13. H reactions with H.^ are an 
essential ingredient of disk chemistry. The decrease in the abundance error bars for Hjj" in the 
model with refined rate uncertainties leads to smaller abundance uncertainties in the related 
species: HCO+, H2CO, H2O, HCN, HNC, HCNH+, and CS. Table [3] contains a number of 
reactions important for the evolution of CO and CO2; e.g., formation of light hydrocarbons, 
neutral-neutral reactions of these hydrocarbons with oxygen, and CO photodissociation. The 
former reactions affect the evolution of C2H, while photodissociation is an important process 
for CN. 

Below we discuss the contributions of each reaction type in detail. 



4-2.1. Radiative association 



The most problematic reactions in the disk are radiative association (RA) reactions. 
Radiative association leads to the formation of a larger molecule from two smaller species 
upon their collision, and requires the emission of excess energy in the form of a photon. 
RA reactions allow the form ation of new bonds and more chemically "advanced" species 



fivan Dishoeck fc Blakd Il998[ ) 



Williams 


1972; 


Herbst 


1980 



Bates fc Herbstlll988l : ISmithI Il989l ). Nonetheless, assigned uncertainties for many of these 
reactions in RATE 06 are only a factor of 2 and less. Since we assume that the rate co- 
efficients of such reactions may vary by an order of magnitude, in accord with the recent 
OSU database, it is not surprising that some RA reactions are identified as among the most 
problematic reactions for disk chemistry. Its importance is also rise due to the fact that we 
start calculation from the atomic initial abundances and use time-integrated correlation co- 
efficients. It means that at early evolutionary time chemistry is dominated by RA reactions 
and, therefore, its time-integrated correlation coefficients increase greatly at early times. 

Two RA reactions are particularly important for the chemical evolution and accumu- 
lation of the abundance uncertainties. First is the reaction between H2 and C, which pro- 
duces CH2 everywhere through the disk at all times. The rate coefficient of this reaction 
is not well known, with with an estimated value of 10"^'' cm^ s~^ at T = 10 — 300 K 



(jPrasad &: Huntress! Il980l : IWoodall et al.l 120071 ). This reaction is relevant for the evolution 
of CO, HCO"'", H2CO, CS, and HCN. Second, the reaction between H2 and C"*" leads to the 
formation of CH^ (important for C"*", HCO"*", C S, NH3, and CN), with an estimated rate 
coefficient of ~ 4 x 10~^^ cm^s~^ (jHerbstlll985l ). These two reactions initiate processes of 
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C and H addition in the chemical model and thus should be relevant for nearly all species 
containing carbon and hydrogen. 



The next most important RA reaction (8th place out of 56) is the formation of pro- 
tonated met hane from molecular hydrogen and CH;^ {a = 1.3 x 10^^^ cm^s^^, (3 = —1; 
SmithI Il989l ). which is an important molecular ion for CO destruction. M ore rece n tly th e 
rate of this reaction has been measured at low temperature of 10 K by iGerlichI (1l995l ). 
1.1 ±0.1 X 10~^^ cm^s~^ for para-H2. 

Much less important (27th and 28th place) are two slow reactions between C and H or 
N producing CH or CN radicals. It is interesting to note that the latter reaction has differen t 
rate coefficients in the OSU and RATE 06 databases JSmith et allbood : Iwoodall et allboO?! ). 
While in osu.2007 thi s reaction rate has no tern perature dependence and no activation barrier 
10~^^ cm^s~^; [Prasad fc Huntress! Il980l ). in RATE 06 the corresponding rate is low, 



[a 

7.9 X 10-^9 



cm" 



s ^, and possesses a small barrier of 96 K (relevant temperature range is 



between 300 and 2 700 K: ISingh fc Andreazzal l2000f ) . Thus the RATE06 and osu.2007 rate 
coefficients for the RA reaction between C and N differ by an order of magnitude in a wide 
temperature range. 

The situation is similar for the RA reaction between molecular hydrogen and CH that 
forms CH3 - an important radical for the evolution of formaldehyde. The corresponding 
OSU.2007 rate at 10 K is 2.5 x 10"^^ cm^ s"^ (a = 3.25 x 10~^^ cm3s-\ (3 = -0.6). In 
contrast, in RATE 06 this reaction has a small barrier of 11.6 K and the 10 K rate value is 
lower by an ord er of magnitude, 1.8 x 10 ~^^ cm^s~^ (a = 5.09 x 10^^^ cm^ s~^, f3 = —0.71, 
and 7 = 11.6 K; iBrownsword et al.l 119971 ). 

In addition to radiative association, we list radiative recombination reactions of C"*" and 
CHjI" in the same group in Table O The radiative recombination of ionized carbon is a key re- 
action for the evolution of C^. Its rate coeffic ient at room temper ature is 4.67 x 10^^^ cm^ s~^ 
with an inverse temperature exponent of 0.6 (IWoodall et al.ll2007l ). The as sumed uncertainty 
of this rate (the factor of 10) in osu.2007 is in contrast to the results of iNahar fc Pradhan 
( 119971 ). who have theoretically predicted this rate value and found that the corresponding 
uncertainty is only ~ 10%. In the RATE 06 database the uncertainty of this rate constant 
is set to 50%. 

Finally, the radiative recombination of CHj^ is the least probable channel in the reactive 
collisions between methyl ions and electrons. The rate of this process is estimated to be 
much higher than for other recombination reactions involving atomic species, particularly in 
cold region s, and might be accurately determined i n laboratory (a = 1.1 x 10"^'' cm'^s^^. 



p = -0.5; Miliar et all Il997l : iLe Teuff et all I2OOOI : IWoodall et al.ll2007f ). In the osu.2007 
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database, this rate has a steeper dependence on temperature, P = —0.7 fjSmith et ahl 12004 ). 

We conclude that RA reactions require a particular attention because their rate un- 
certainties are often large and not well known so that the modeling results can be heavily 
affected. 



4-2.2. Cosmic ray ionization 



The next most important reactions that contribute strongly to the resulting abundance 
scatter of many molecules are the cosmic ray ionization of molecular hydrogen and helium 
atoms. The cosmic ray particles able to penetrate deeply into the disk are mostly high-energy 
protons (E >, 100 Mev) and heavy nuclei (e.g. iron, E' ~ 1 Gev) (jPalgarno fc McCraylll972l : 



Leger et al.l Il985l : iDolginov fc Stepinskil Il994j ). Direct ionization initiate chemical evolu- 
tion in dark and cold regio ns where cosmic ray particles remain the main ionization source 
( iHerbst fc Klempererlll973l ). In the disk this is the outer midplane at r > 30 AU and the 
lower part of the intermediate layer (jGammidll996l : ISemenov et al.ll2004l ). Moreover, since 
we added the X-ray ionization rate to the CRP rate value in an attempt to crudely treat the 
effect of high-energy stellar radiation, it is not surprising that high correlation coefficients 
for some CRP-ionization reactions are reached in the disk atmosphere and intermediate 
layer at r < 100 AU for several key species such as CO (see Fig. [^. A more accurate 



approach to model X-ray che r nistry in disks is req uired ( Neufeld et al.l 1 1994 : iMaloney et al. 
19961 : iLepp fc Dalgarnolll996l : Istauber et aPboOsh . 



It is natural that because the cosmic ray ionization of H2 and He is the primal chemical 
process, it is also the key factor for the accumulation of abundance scatters. In the RATE 06 
da tabase the correspon ding rate coefficients are estimated to be accurate within a factor of 
2 (jWoodall et al.l 120071 ). However, this rate uncertainty has two origins. 



First, there is the physical uncertainty that is related to our limited knowledge of the 
absolute CRP flux and its penetrating ability, which may vary in different a strophysical 
objects ( ISpitzer fc Tomaskdll968l : lDolginov fc Stepinskilll994j : ICaselli et al.lll998l ). While the 
direct ionization process driven by CRPs is well stu died and the corresponding cross sections 
and rates are well known (ISolomon fc Wernerlll97ll ). the second chemic al uncertainty resides 
in th e treatment of the ionization by energetic secondary electrons (IGlassgold fc Langei 
19731 ). Here the elemental composition of the gas plays a decisive role (jllgner fc NelsonI 



200 



i. 



The CRP-induced ionization rate of H2 currently used in all UMIST and OSU databases 
was taken from the PhD thesis of Black (1975). It was computed for a mixture of atomic 
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and molecular hydrogen but new calculations for a mixture of hydrogen an d helium would 



be required in order to have a more consistent rate value (see discussion in Wakelam et al. 



2006 



Other CRP-ionization reactions listed in Table [3] do not have a reference in the 
RATE 06 database so it is hard to judge how accurate their rates are. 

Apparently, the rate uncertainty of the CRP-driven reactions can be partially eliminated 
by more accurate modeling of the full cascade of ionizing event s, but the physic a l uncer tainty 
will still remain, which needs more efforts from observers. I Wakelam et al.l (j2006bl ) have 
demonstrated how crucial this may be for the results of pure gas-phase chemical models of 
static molecular clouds. 



4.2.3. Photoreactions 



Another set of chemical processes with problematic rates are photodissociation (PD) and 
photoionization (PI) reactions, as listed in Table [31 These reactions are important in those 
regions where UV photons are either able to penetrate (the disk atmosphere and inter mediate 
layer) or are induced by cosmic ray particles (disk midplane: [Prasad fc Tarafdarlll983l ). High- 
energy UV photons destroy or ionize gas-phase molecules and allow surface species to desorb 
back into the gas phase. 

While both the PD and PI cross sectio ns for many species can be measured or calcu - 
lated with a rather good accuracy of < 50% (Ivan Dishoecklll988l : Ivan Dishoeck et al.ll2006l ). 
the major difficulty in obtaining accurate photorates is our limited knowledge of the UV 
radiation field i nside protoplanetary disks . 

citetBeaOS and Ivan Dishoeck et al.l (120061 ) have shown that many T Tauri stars emit strong, 
non-blackbody UV radiation. Their spectra resemble that of the interst ellar (IS) UV radia- 



tion fi eld, with a large fraction of the UV fiux emitted in the Lya line. Ivan Zadelhoff et al. 



( 120031 ) have demonstrated the importance of scattering for the UV penetration deep into 
the disk interior. Furthermore, for the two most abundant molecules, H2 and CO, that ef- 
fectively dissociate through partly overlapping lines shortward of about 1 120 A, the optical 
thickness in t hese lines can become so high that self- and mutual-shielding have to be taken 



into account (IDalgarno fc Stephenslll970l : lvan Dishoecklll988l : lLee et al.lll996l : IJonkheid et al. 
2006h . 



In the U MIST and OSU datab ases, the photorates have m ostly been adopted from the 
compilation of Ivan Dishoeck! (Il988l ) and iRoberge et al.l (Il99ll ) , where they have bee n com- 
puted using a simple plane-parallel approach and the IS (IS) UV radiation field of iDraine 
( 119781 ). The IS UV fiux is characteristic of a diluted radiation field from an early B star 
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(Spitzer 


1978; 


1 — , ^^^^ 

van Dishoeck 


1994) 



resentative of the UV radiatio n from much cooler Herbig Ae and T Tauri stars surrounded by 



circumstellar disks. Recently, Ivan Dishoeck et al.l (120061 ) have recalculated the correspond- 



ing PD and PI rates, using more appropriate UV spectral distributions typical of Herbig Ae, 
T Tauri, and Sun-like stars. 

The 10 problematic photodissociation and photoionization reactions found by the sen- 
sitivity analysis are involved in either the formation of key species (C"*" and CN), their 
destruction (CO, CN, CS, NH3, and HCN), or the destruction of parental molecules (OH, 
CH, and CH2). Their rates in the RATE 06 database are in general inaccurate by a factor of 
2 with some exceptions. The photodissociation of CO proceeds indirectly, through predis- 
sociation, and depends on the amount of molecular hydro gen and CO molecules in the line 
of sight to the source of UV radiation (Ivan Dishoeckl 119871 ). and thus has a rate uncertainty 
of a factor of 10 (see Table [3]). In contrast, the rates of direct photodissoc iation of NH3, 
OH, and HCN are inaccurate by the factor of 1.5 only (IRoberge et al.lll99ll ). Despite such 
small error bars, the latter three reactions were identified as problematic reactions for disk 
chemistry. 



4-2. 4- Ion-neutral and neutral-neutral reactions 



Ion-neutral (IN) and neutral-neutral (NN) reactions form the largest fraction of the 
problematic reactions (32 out of 56) for disk chemistry, but their individual importance is 
not as significant as for CRI and RA reactions (Table [3]). 

Ion-neutral reactions are usually exothermic and rapid, and lead to the formation of 
new species by bond rearrangement. Their rate coefficients can often be obtained by simple 
Lange yin theory, which relates the isotropic polarizability of the species and its reduced 
mass (IClarylll988l ). For molecules that possess a large dipole moment, the rate values can be 
significantly enhanced at low temperature by the long-dista nce Coulomb attraction betwee n 
the positive ion and negatively charged side of the molecule ( lAdams et al.lll985l : IClarylll985l ). 
Moreover, many neutral-neutral reactions with small barriers or without a barrier can also be 
fast under interstellar conditions. The rates of most radical-radical reactions and even some 
radical-stable neutral reactions are contro lled by long-range attractive forces and thus do not 



get sm aller at low temperatures (see, e.g.. lSmithlll988l : ISims et al.lll993l . ll994bl : ICanosa et al. 



19971 ). The typical rate coefficient for these reactions is about an order of magnitude lower 
than for ion-neutral processes. The rates of some of the fast neutral-neutral reactions at 
low temperatures have been measured in the laboratory and predicted theoretically (e.g.. 
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Clary et al.lll994J : ISims et al.lll994al : IChastaing et al.ll200ll : ISmith et alJl2004l ). 



In RATE 06, most of the IN and NN rate coefficients are uncertain by a factor of 2 
and less, up to the accuracy of 25% (Table [3]). However, a major ambiguity that resides 
in the IN rate constants is the treatment of their temperature dependence. In two versions 
of the RATE 06 database the rate value for the reaction involving a polar molecule and 
a molecular ion can have a negative dependence on temperature {(3 = —0.5; adopted in 
this work) or no temperature dependence {(3 = 0), while its rate coefficient fits the value 
at room temperature. In contrast, in the OSU network the corresponding IN rates based 
on low temperature theoretical estimations to be most accurate at ~ 10 K even though 
the listed reaction rate coefficients are scaled with respect to 300 K with /5 ~ —0.5. Two 
problematic reactions involving a polar molecule and an ion with significantly different rates 
in RATE 06 and osu.2007 are those between C"*" and NH, forming CN+ and H (aosu.2007 = 
4.6 X 10~^ cm^s~^ vs. ctRATEoe = 7.8 x 10~^° cm^s~^) and between and CH, producing 
CS+ and H fansn 2007 = 4.4 x 10"^ cm^s"^ vs aRATRnfi = 6.2 x 10"-^° 



cm'^ s ^ ] 



Other IN reactions with different rates in the RATE 06 and osu.2007 databases are 
reactions involving an atom and a molecular ion. The break up of molecular hydrogen upon 
collision with helium ion involves a small barrier of 35 K and a = 3.7 x 10~^^ cm^s~^ in 



RATE 06 (the relevant temperature range is between lOK and 300K; IWoodall et al.l 120071 ). 
while in osu.2007 the same rate coefficient is about 5 times smaller and has no barrier 



(jSmith et al.ll2004l ). However, this results in difference of the osu.2007 and RATE 06 rate 
values of an order of magnitude at most. 



Similarly, the OSU reaction between molecular hydrogen and ionized ammonia produces 
hydrogen atom and protonated ammonia with a rate coefficient of 1.5 x 10~^^ cm^s~^ and 
a steep temperature dependence (/3 = —1.5). In RATE 06 this reaction has a non-Langevin 
rate with /? = 0. At temperatures below 20 K its rate coefficient i s twice as large, a = 3.36 x 



10 cm s . and a negative activation barrier of -35.7 K exists (lAdams fc Smitlj|l984j ). At 



temperatures between 20 and 300 K this reaction has no barrier and a = 2 x 10~^^ cm^s~^. 
For temperatu res above 300 K the r ate coefficient is 1.7 x 10^^^ cm^ s^^ and a barrier is 
about 1 000 K (IFehsenfeld et al.lll975l ). Consequently, the osu.2007 and RATE 06 rate values 
can differ by up to one order of magnitude. 

Finally, the IN reaction between H.^ and O produces only 0H+ and H2 in OSU, but in 
RATE 06 a second, less probable alte rnative channel that leads to the production of ionized 
water and atomic hydrogen is given (IMilligan fc McEwaru I2OOOI ) . 



Note that the 4 IN reactions with the smallest rate uncertainty of 25% are still prob- 
lematic for disk chemistry. Such high sensitivity of the disk modeling results to the rates 
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of these reactions is caused by their importance for the evolution of CO, HCO^, and HCN. 
The most important one is the IN reaction between hehum ions and CO, which leads to 
the hypersensitivity of the final abundances of many carbon-bearing species in the inner 
part of the disk intermediate layer (see Fig. [2]). The other reactions include the pri- 
mal route to the formation of HCO"*" and destruction of CO: Ht -|- CO HCO"*" -|- H2 



Q S — 1 

cm s ■ 



[a = 1.7 X 10 
by ionized helium {a 
by ionized carbon {a ■ 



Kim et al.lll975l ) as well as the destruct ion of molecular nitrogen 



6.4 X 10"^° cm^s"^; IWoodall et aD l2007h and destruction of HCN 



3.1 X 10^9 cm3s~\ 



-0.5; 



Clarv et al.lll985h . 



All of the 15 problematic neutral-neutral reactions for disk chemistry are reactions 
involving atoms (O, N, C, H, S) and either light hydrocarbons (CnHm, n=l,2, m=l-3) or CN, 
HCO, and H2CO. In contrast to the IN reactions, the NN reactions are mainly important for 
the chemical ev olution of the key species in warm disk regions, in particular reaction with 
atomic oxygen (ITscharnuter fc Gail! 120071 ). Though their rates are rather accurate, with 
uncertainties that are typically not higher than 50%, some of the NN rates are still different 
in RATE 06 and osu.2007 databases. 

The NN reaction between N and HCO produ ces only OCN and H in osu.2007, with a 
rate coefficient of 10~^° cm^ s~^ (jSmith et al.ll2004l ). In RATE 06 this reaction has in addition 



2 other channels that form either HCN and O (a = 1.7 x 10 ^° cm^s ^) or CO and NH 



(a 



5.7 X 10-12 



cm^ s 



\ (3 = 0.5, and 7 = 1 000 K; IWoodall et al.ll2007f ) 



Another reaction in the list is H + CH, which leads to C and H2. In osu.2007, the rate 
coefficient has a weak depende nce on temperature and no barrier {a = 2.7 x 10-^^ cm'^s^^, 
f3 = 0.38; lie Teuff et aP boOoh. while in RAT E 06 this rate has /3 = and an 80 K barrier 
(a = 1.3 X 10 



cm^ s 1 ; 



Woodall et all 120071 ). 



Finally, the neutral- neutral reaction between H and CH2, which forms CH and H2, has 

It c 

this value is only 6.6 x 10^^^ 



the OSU.2007 rate coefficient of 2.7 x 10 cm^ s ^ taken from the UMIST95 database, but 



cm'^ s 



in RATE 06. 



We conclude that the reaction rates of many important ion-neutral and neutral-neutral 
reactions are not known accurately enough, as a major controversy in their temperature 
dependence still persists. This is particularly true for the reactions involving an ion and 
a polar molecule, whose rates in RATE 06 can differ by a factor of ~ 5 at 10 K for the 
dipole/non-dipole versions of this ratefile. The ion-polar rates in osu.2007 can be larger 
than those in the dipole version of RATE06 by another factor of 4-5, especially for linear 
neutral reactants, in which allowance is made for sub-thermal rotational populations via the 
so-called "locked- dipole" approach. 
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4-2.5. Dissociative recombination 



Molecular ions are efficiently converted into other, less complex molecules via disso- 
ciative recombination (DR) with electrons or charged grains. Similarly to IN reactions, 
these processes are especially fast at low temperatures with typical rate coefficients of about 
10"^ cm~^ s~^ at 10 K due to the long-range Coulomb attraction. The DR rates can be mea- 
sured rather accurately with afterglow and storage ring techniques ( iFlorescu-Mitchell fc Mitchell 
20061 ) . However, the sensitivity of their branching ratios with respect to the tempera- 
ture is mostly unexplo red and may not be described by a simple power law at low T 
(jPetrignani et al.ll2005l ). 



A strong influence of the DR rat es and branching ratios on the results of chemical models 



of mol ecular clouds has been found (IMillar et al.lll988l : ISemaniak et al.ll200ll : iGeppert et al. 
2005al ). Our sensitivity method picked up two DR reactions that introduce significant un- 
certainties in computed abundances in protoplanetary disks. These include dissociative 
recombination of protonated ammonia into ammonia and hydrogen atoms (50th place out 
of 56), and a main destruction channel for HCO"^ (19th place, see Table [3]). Note that their 
rates are supposed to be known within an uncertainty of 25%. 

The DR of HCO^ in RATE 06 has one dissociative channel to CO and H with the total 



rate coefficient of 2.4x lO"'^ cm^ s'^ and f3 = -0.69 (lMitchelllll990r i. In the osu.2007 database 
2 additional channels are listed, which lead to the formation of either OH and C or CH and 
O but with a much smaller probability of ~ 4% each. The DR of NH^ in both RATE 06 
and OSU.2007 involves 3 branching channels into NH2 and either two hydrogen atoms (21%) 
or one molecular hydrogen (10%), or am monia and atomi c hydrogen (69%) with the total 
rate of 1.5 x lO"*' cm^ s"^ and (3 = -0.5 Jvikor et al.lll999h . 



Still, more effort needs to be invested into the (re-)investigation of the DR branching 
ratios and their temperature-dependence, especially at low ~ 10-20 K temperatures. 



5. Summary and conclusions 

The influence of gas-phase reaction rate uncertainties on the results of disk chemical 
modeling has been studied by a Monte Carlo method. The rate coefficients in the RATE 06 
network were varied 8 000 times within their uncertainty limits using a log-normal distribu- 
tion. Sets of abundances, column densities, and their error bars were computed for a number 
of key species in protoplanetary disks. We found that typical uncertainties of the molecular 
column densities do not exceed a factor of 3-4 even for the largest of the key molecules, 
which is comparable with observational uncertainties. The column densities of CO, C"*", 
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, H2O, NH3, N2H+, and HCNH+ have particularly small error bars. A straightforward 
correlation analysis between molecular abundances and reaction rates was performed for the 
entire evolutionary time, and the most problematic reactions involving CO, C^, CS, HCO"'", 
H2CO, CN, HCN, and NH3 were identified. We showed that the rate coefficients of about a 
hundred chemical reactions constituting only a few percent of the entire RATE 06 database 
need to be determined more accurately in order to significantly decrease uncertainties in 
the modeled abundances and column densities of the key observable molecules. We argue 
that it is worthwhile to (re-)investigate, either experimentally or theoretically, the rate co- 
efficients of basic radiative association and cosmic ray ionization reactions as well as the 
temperature-dependence of key ion-neutral and neutral-neutral reactions, and the branching 
ratios and products of dissociative recombination processes. The rate uncertainties of cosmic 
ray ionization and photoreactions are partly due to uncertain physical parameters such as 
the CRP flux and UV penetration, and thus can be decreased with better physical models 
and observational data. 
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R, AU R, AU 

Fig. 1. — Temperature (K) and density (cm~^) structure of the adopted disk model. The 
total scale height of the disk Zmax increases with the radius such that at r = 100 AU it has 
a value of 90 AU, whereas Z^^ « 1 300 AU at r = 800 AU. 
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Fig. 2. — Histograms of the abundance distributions show the number of chemical models in 
which a given molecular abundance is reached at 5 Myr. (Top) Non-Gaussian, double-peak 
distributions are shown for HCO"*" (left) and CO (right) in the inner disk region (r ^ 97.5 AU, 
z 30 AU). (Bottom) Nearly Gaussian distributions are shown for the same species but 
located farther out from the central star (r 382 AU, z 205 AU). The point numbers in 
the boxes refer to the disk cells for which these calculations were performed. 
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R, AU 

Fig. 3. — Distribution of mean molecular abundances relative to the total amount of hydrogen 
nuclei for several observationally important species in the disk at 5 Myr. The distributions 
are obtained by averaging the 8 000 model results calculated with randomly varied reaction 
rates. Results are shown for CO, C"*", CS, HCO"*" (top row, from left to right) and H2CO, 
NH3, HCN, and CN (bottom row) at r > 50 AU. The disk vertical axis is given in units of 
the pressure scale height. 
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Fig. 5. — The same as in Fig. H] but for radial distributions of the column densities. The 
central solid line represents the model-averaged column densities. The corresponding column 
density scatters computed with the standard model are depicted with dotted lines. The 
column density scatters calculated with the chemical model where the rates of the 56 most 
problematic reactions are varied by small factors of 1.25-2 only are shown with solid lines. 
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Fig. 6. — The largest correlation coefficients of relevant reactions for the evolution of NH3 
as a function of time in the disk intermediate layer (r = 97 AU, 2; = 30 AU). 
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Table 1. Initial abundances 



Species 


n(X)/n(H) 


He 


9.75(-2)^ 


H2 


4.99(-l) 


H 


2.00(-3) 


C 


7.86(-5) 


N 


2.47(-5) 





1.80(-4) 


s 


9.14(-8) 


Si 


9.74(-9) 


Na 


2.25(-9) 


Mg 


1.09(-8) 


Fe 


2.74(-9) 


P 


2.16(-10) 


CI 


1.00(-9) 


F 


I.OO(-IO) 


^A(-B) 


means 



AxlO-^ 
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Table 2. Typical column density uncertainties 



Species 


Standard 


Improved 




model 


model 




factor ~2 


~1.5 


c+ 


factor ~2 


~1.5 


CO 


factor ~1.25 


1.1 


C02 


factor ~3.5 


~2.5 


cs 


factor 3 


2.5 


HCO+ 


factor ~ 3 


~1.5 


H2CO 


factor 3 


~1.5 


CN 


factor 4 


2 


HON 


factor ~ 3 


2.5 


HNC 


factor 3 


~ 2.5 


N2H+ 


factor 2.5 


2 


NH3 


factor 2.5 


2.5 


HCNH+ 


factor 2.5 


2.5 


H2O 


factor 1.7 


1.4 


C2H 


factor 4 


~ 3 
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Table 3. Most problematic reactions 



Reaction 



Uncertainty Weight Type 



factor 10 


1.0 




RA 


factor 10 


4.2( 


-1) 


RA 


factor 10 


1.6( 


-2) 


RA 


factor 10 


7.6( 


-2) 


RA 


factor 10 


2.3( 


-2) 


RA 


factor 10 


2.2( 


-2) 


RA 


factor 10 


1.2( 


-2) 


RR 


factor 10 


1.2( 


-2) 


RR 


factor 2* 


1.4( 


-1) 


CRI 


factor 2 


1.0( 


-2) 


CRI 


factor 2 


2.3( 


-1) 


CRI 


factor 2 


9.7( 


-2) 


CRI 


factor 2 


3.8( 


-1) 


PI 


factor 2 


6.7( 


-2) 


PI 


factor 10 


5.6( 


-2) 


PD 


factor 2 


2.3( 


-2) 


PD 


factor 2 


6.8( 


-3) 


PD 


factor 2 


3.2( 


-2) 


PD 


factor 2 


2.6( 


-2) 


PD 


factor 1.5 


8.2( 


-3) 


PD 


factor 1.5 


1.5( 


-2) 


PD 


factor 1.5 


7.8( 


-3) 


PD 


factor 1.25 


9.6( 


-3) 


IN 


factor 1.25 


7.4( 


-3) 


IN 


factor 2 


3.3( 


-2) 


IN 


factor 1.5 


1.0( 


-2) 


IN 


factor 1.25 


7.3( 


-3) 


IN 


factor 2 


4.3( 


-2) 


IN 


facjtor 2 


2.0( 


-2) 


IN 


factor 2 


2.1( 


-2) 


IN 


factor 1.5 


4.4( 


-2) 


IN 


factor 1.5 


8.4( 


-3) 


CT 


factor 2 


1.3( 


-2) 


IN 


factor 2 


3.1( 


-2) 


IN 



H2 + C 
H2 + C+ 

Ha + CH 
H2 + CH+ 
C + H 

C + N 
C+ + e- 
CH+ + e- 

H2 + C.R.P. 

H2 + C.R.P. 
He + C.R.P. 
N + C.R.P. 

C + hi^ 

C + C.R.P./11/ 

CO + hiy 

CH + hiy 
CB.2 + hv 
CN + hu 
CS + hv 
NH3 + hv 
HCN + hv 
OH + hv 

He+ + CO 
He+ + N2 
H3+ + C 
O 

CO 
H+ +N 

H2 + He+ 
H2 + He+ 
H2 + NH+ 
H+ + O 
C+ + CH 
C+ + OH 



CH2 + \iv 
CHj + \iv 

CH3 + hv 
CH+ + hv 
CH + hv 
CN + hv 
C + hv 
CH3 + hv 

H+ + e- 

H+ + H + e- 
He+ + e- 
N+ + e- 

C+ + e- 
C+ + e- 
O + C 
C + H 
CH + H 
N + C 
S + C 
NH2 + H 
CN + H 
O + H 

O + C+ + He 
N + N+ + He 
CH+ + H2 

0H+ + H: 
HCO+ 
NH+ + H 

He + H+ - 
He + H+ 
NH+ + H 
0+ + H 
C+ + H 
C0+ + H 



I2 
H2 



H 



Table 3 — Continued 



Reaction 



Uncertainty Weight Type 



C+ + NH 
C+ + HCN 
C + HCO+ 
CH + S+ 
CH+ + O 
NH^ + Mg 

C + CH2 
C + C2H 
N + CH2 
N + CN 
N + HCO 
H + CH2 
H + CH 
CH + O 
CH + O 
CH + N 
CH + S 
CH2 + S 
CH2 + o 
CH3 + o 

HCO+ + e- 

Nn+ + 0- 



CiN+ H 




l.li 




IN 


CNC+ + H 


factor 1 25 


7.9( 




IN 


CO 4- CH+ 


Tfl oi" nv 






IN 


r;s+ -1- H 




Z.\Ji 


-21 


IN 


HCO+ -1- H 


idUUJi Zi 


7 9^ 




TN 


NH3 + Mg+ 


factor 2 


i.i( 


-2) 


CT 






3.6( 


-21 


NN 






0.0 


•^l 


NN 


HCN + H 


factor 1.5 


8.5( 


-2) 


NN 


N2 + C 


factor 2 


4.6( 


-2) 


NN 


HCN + 


factor 2 


2.9( 


-2) 


NN 


CH + H2 


factor 1.25 


1.9( 


-2) 


NN 


C + H2 


factor 1.5 


1.1( 


-2) 


NN 


HCO+ + Q- 


factor 1.5 


4.0( 


-2) 


NN 


CO + H 


factor 1.25 


8.3( 


-3) 


NN 


CN + H 


factor 1.25 


2.3( 


-2) 


NN 


CS + H 


factor 2 


1.6( 


-2) 


NN 


CS + H2 


factor 2 


6.3( 


-2) 


NN 


CO + H + H 


factor 1.25 


1.0( 


-2) 


NN 


H2CO + H 


factor 1.25 


1.4( 


-2) 


NN 


CO + H 


factor 1.25 


3.1( 


-2) 


DR 


Nil:; + 11 


fad or 1.2") 


<S.1( 




DR 



* These rates are partly uncertain due to uncertainties in physical 
parameters. 



Note. — (CRI) Cosmic ray ionization; (CT) Charge transfer; (DR) 
Dissociative recombination; (IN) ion-neutral reaction; (NN) Neutral- 
neutral reaction; (PD) Photodissociation; (PI) Photoionization; (RA) 
Radiative association; (RR) Radiative recombination. Three most im- 
portant reactions for the disk chemistry as identified by the sensitivity 
analysis are written in boldface. 



